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Models for gamma-ray production in low-mass microquasars 
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Unlike high-mass gamma-ray binaries, low-mass microquasars lack external sources of radiation and matter that could 
produce high-energy emission through interactions with relativistic particles. In this work we consider the synchrotron 
emission of protons and leptons that populate the jet of a low-mass microquasar. In our model photohadronic and 
inverse Compton (IC) interactions with synchrotron photons produced by both protons and leptons result in a high- 
energy tail of the spectrum. We also estimate the contribution from secondary pairs injected through photopair 
production. The high-energy emission is dominated by radiation of hadronic origin, so we can call these objects 
'proton microquasars'. 

Keywords: X-rays:binaries - gamma rays: theory - radiation mechanisms: non-thermal 



1. Introduction 

Microquasars are accreting binary systems that 
accelerate particles up to relativistic energies, as 
demonstrated by the detection of synchrotron radi- 
ation from radio up to X-ray energies. In micro- 
quasars with a high- mass donor star, the relativis- 
tic jets ejected from the surroundings of the com- 
pact object must traverse photon and matter fields 
produced by the star. This results in interactions 
that can produce high-energy gamma rays, as dis- 
cussed by several authors.'^'"' In fact, three high- 
mass microquasars have already been detected as 
gamma ray emitters. Low-mass microquasars, on 
the contrary, have not been detected so far. High- 
energy emission from these objects has been mod- 
eled in Refs. 4, 8, 9 and 10. These authors have fo- 
cused only on leptonic processes. Here we present 
a new family of models for the non-thermal emis- 
sion of low-mass microquasars which include several 
improvements. Among other things, we consider an 
hadronic content in the jets, proton synchrotron ra- 
diation, photomeson and photopair production, lep- 
tonic processes for primary and secondary particles, 
etc. Relativistic particle distributions are calculated 
in self-consistent way, solving the transport equation 



under reasonable assumptions. Our models predict 
that low-mass microquasars might be detectable by 
instruments like GLAST and Chcrcnkov telescopes 
like VERITAS, HESS and MAGIC, and other future 
extensions of these arrays. 

2. Model 

We consider an inhomogeneous jet of conical geome- 
try perpendicular to the orbital plane of the binary. 
The jet is launched at a distance zg from the com- 
pact object and has an initial radius — Q.l zq. The 
outflow carries a significant fraction of the accretion 
power, ^^'^^ Ljct — gjot-^accr- We further assume that 
a fraction = ^Ljot is in the form of relativistic 
protons and leptons, Lroi — Lp + Le- We relate the 
energy budget of both species as Lp — aL^. Parti- 
cles are accelerated by diffusive shock acceleration 
leading to a power-law injection 



erg ^s ^cm ^. 



(1) 



We consider a canonical value a = 2.2 for the spec- 
tral index, and we assume that the injection rate 
decreases with the distance to the compact object. 
Acceleration of particles takes place in a compact re- 
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gion centered at z^cc^ with a thickness Az = 5r (-Zacc)- 

The cfScicncy of the acceleration mechanism is char- 
acterized by the parameter rj, in such a way that the 
acceleration rate is 



Tabic 1. Values of the various parameters cliaractcriziiig the jet and 
the distributions of relativistic particles that were used in the calcula- 
tions. 



*acc = VceBE- 



(2) 



We assume an efficient accelerator with t] = 0.1. The 
magnetic field in the jet decreases as the jet expands 
as B = Bq{zq/z). Its value at zq is determined by 
equipartition between magnetic and kinetic energy 
densities, Ub{zo) = Ul^^{zo). This yields Bq ~ fO^ G 
and makes synchrotron losses the dominant channel 
of radiative particle cooling. For a particle of mass m 
and energy E, the synchrotron cooling rate is given 

by 



Parameter 


Value 




LuccT- accretion power 


1.7 X 10^9 erg 


s-iW 


Qjct- coupling constant accretion-jet power 


0.1 




zo: jet's launching point 


2 X 10*^ cm(2) 




Zacc- location of the acceleration region 


10^ cm 




Az: thickness of the acceleration region 


5r(Zacc) 




Fjct ; jet's bulk Lorentz factor 


1.5(1) 




0: viewing angle 


30° 




4: ratio L^^l/Ljet 


0.1, 10-3 




a: ratio Lp/Le 


1, 100, 103 




r): acceleration efficiency 


0.1 




-B(): magnetic field at za 


2 X 10^ G 




i?™g": minimum particle energy 


100 mp,eC^ 




gmax. maximum particle energy p,e~ 


lOio, 7 X 1016 


eV(3) 



Typical value for the low-mass microquasar XTE J1118+480l* 
50Rschw for a 8M0 black hole as XTE J1118-|-480i5 
(3) These are the maximum values obtained along the jet. 



t. 



synchr 



4 / me 
3 V m 



I 3 cotUb E 



(3) 



Particle distributions in the steady state 
Np,e{Ep,e,z) (erg-^cm-"^) are calculated solving a 
transport equation that takes into account parti- 
cle injection, cooling and particle escape from the 
acceleration region, 



3. Results 

In addition to synchrotron radiation, we consider 
other two photon production mechanisms: IC and 
proton-photon colissions (pj) on the synchrotron ra- 
diation fields of both protons and leptons. The later 
process has two main branches: photomeson produc- 
tion 

J) + 7 — > J) + a-K^ + b (tv~^ + 7r~) 



dE 



dE 
'dt 



N{E,z) 



N{E,z) 



Q{E,z), (4) 



where dE/dt\-^^^^ is the sum of all radiative and me- 
chanical losses for each type of particle. 
We estimate the escape time as tesc = ^acc/^jet) 
where fjet ~ 0.7 c corresponds to a bulk Lorentz fac- 
tor Fjet = 1.5, typical of a jet in the low-hard state 
of the micro quasars. The proton spectrum goes as 
Np oc S"^-^ for all Ep. In the case of electrons, in- 
stead, synchrotron losses introduce a steepening that 
leads to N^. oc E~^ '^ at all energies before the cutoff. 
The acceleration and cooling rates for both protons 
and electrons at the acceleration region are shown in 
Figure 1 . The maximum energy of all particles is de- 
termined by the synchrotron cooling rate. In Table 
1 we show some representative values of the relevant 
parameters. 



and photopair production 

The spectrum from the decay of 7r° 2j is esti- 
mated following Ref. 16. For the injection of elec- 
tron/positron pairs we use the formulae given by 
Refs. 17 and 18. The IC spectra are calculated in 
the local approximation of Ref. 19, whereas for syn- 
chrotron radiation we use the classical expressions, 
as given, for example, by Ref. 20. All the calcula- 
tions are performed in the jet co-moving reference 
frame and then transformed to the observer frame 
using the appropriate Doppler factor. We &x 6 = 30° 
as the viewing angle for numerical estimates. Since 
the jet is only mildly relativistic, none of our results 
depends strongly on this value as long as it is not 
very close to zero. 
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Fig. 1. Acceleration and cooling rates at Zacc for electrons (left panel) and protons (right panel), for a = 1. In the case of protons, 
the p7 cooling rate is the sum of the photopair and photomeson production contributions. The dominant channel of energy loss 
for electrons is synchrotron radiation. For protons, adiabatic (low energies) and synchrotron (high energies) losses are the main 
cooling mechanisms. 



Figure 2 (b) , (c) and (d) correspond to a proton- 
dominated jet. The peak of the spectral energy dis- 
tribution in the three cases is due to the proton syn- 
chrotron radiation. At energies beyond 10 TeV, the 
dominant contribution comes from the decay of neu- 
tral pions produced through photomeson mechanism. 
In the energy range between 1 TeV and 100 TeV, the 
synchrotron radiation from photopairs is systemati- 
cally more significant than that generated by sec- 
ondary pairs from charged pion decays. 

Case (a) corresponds to equipartition between 
protons and leptons (a = 1). As expected, the lepton 
synchrotron luminosity rises. A denser photon field 
results in an enhanced (although relatively small) 
leptonic inverse Compton (IC) emission, but an ap- 
preciable P7 luminosity at very high energies, charac- 
terized by a hard spectrum. In case (d) the relativis- 
tic particle content is smaller, ^ ~ Lioi/ijct = 10~^. 
This reduces synchrotron, IC and P7 luminosities by 
one and two orders of magnitude, respectively. 

Concerning the radio emission, our models, 
which have a high minimum injection energy for 
electrons, are not intended to reproduce the observ- 
able radio spectrum of systems like XTE J1118+480 
(^radio ~ iQ^^"^^ crg s^^^^). Howcvcr, modcls with 
iJ™'" ~ 2meC^ can yield the expected luminosity and 
spectrum in this band.^^ 



4. Conclusions 

Low-mass microquasars with steady jets can be sig- 
nificant sources of high-energy radiation if they have 
an important hadronic content. These objects can be 
found out of the galactic plane, because they have 
large proper motions. They can also be sources of 
neutrinos through charged pion decays from p7 inter- 
actions, as already noticed by Ref. 24. These 'proton 
microquasars' might be detected in the near future 
by GLAST and ground-based Cherenkov telescope 
arrays, since gamma-ray luminosities in the range 
1033 _ 1036 erg g-i between 1 GeV and 10 TeV can 
be expected under a variety of conditions. Moreover, 
some of the variable, unidentified sources detected 
by EGRET around the galactic center at > 100 
MeV could be 'proton microquasars' ejected from the 
halo or the plane of the Galaxy. 
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Fig. 2. Production spectra obtained for different values of jet content of relativistic particles 5, and proton-to-lepton luminosity 
ratio a. 



References 

1. Mirabel I.F. et al, Nature 358, 215 (1992) 

2. Corbel S. et al, Science 289, 196 (2002) 

3. Romero G.E. et al, A&A 410, LI (2003) 

4. Bosch- Ramon et al, A&A 447, 263 (2006) 

5. Aharonian F.A. et al., Science 304, 746 (2005) 

6. Albert J. et al. Science 312, 1771 (2006) 

7. Albert J. et al, ApJ 665, L51 (2007) 

8. Atoyan A.M. & Aharonian F.A., MNRAS 302, 253 
(1999) 

9. Markoff S. et al, A&A 372, L25 (2001) 

10. Grenier LA. et al, Ap&SS 297, 109 (2005) 

11. Falcke H. & Biermann P., A&A 293, 665 (1995) 

12. Kording E.G., Fender R.P., & Migliari S., MNRAS 
369, 1451 (2006) 

13. Fender R.P., Belloni T.M., Gallo E., MNRAS 35, 



1105 (2004) 

14. Markoff S., Falcke H., & Fender R., BAAS 32, 1546 
(2000) 

15. Gelino D.M. et al, ApJ 642, 438 (2006) 

16. Atoyan A.M. & Dermer C. D., ApJ 586, 79 (2003) 

17. Chodorowski M.J. et al, ApJ 400, 181 (1992) 

18. Mastichiadis A. et al, A&A 433, 765 (2005) 

19. Ghisellini G., Maraschi L., & Treves A., A&A 146, 
204 (1985) 

20. Blumenthal G.R. & Gould R.J., Rev. Mod. Phys. 42, 
237 (1970) 

21. Fender R.P., Hjellming R.M., Tilanus R.P.J, et al, 
MNRAS 322, L23 (2001) 

22. Romero G.E. & Vila G.S., A&A, submitted (2008) 

23. Mirabel I.F. et al. Nature 413, 139 (2001) 

24. Levinson A. & Waxman E., Phys. Rev. Lett. 87, 
171101 (2001) 



